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SUMMARY 


This report contains the results of a study developed for the Langley 
TCV B-TS? program to apply existing navigation aids plus magnetic leader 
cable signals and develop braking and reverse thrust guidance laws to pro- 
vide for rapid automated rollout, turnoff, and taxi to reduce runway occu- 
pation time for a wide variety of landing conditions for conventional 
commercial-type aircraft. 

The study concludes that it is possible and feasible to execute rollout 
and turnoffs on conventional runways in under 25 seconds with an integrated 
automated landing deceleration system for a large variety of aircraft and 
runway landing conditions. 

Closed loop guidance laws for braking and reverse thrust are derived 
for rollout, turnoff, and taxi, as functions of the landing speed, the desired 
taxi speed and the distance to go. Brake limitations for wet lunway conditions 
and reverse thrust limitations are taken into account to provide decision rules 
to avoid tire skid and to choose an alternate turnoff point, farther down the 
runway, to accommodate extreme landing conditions. 

The ground rules underlying the guidance policy are: 

1) To exit the rollout runway in under 25 seconds, if possible; 

2) The guidance shall be closed loop, using the brakes as a modulating 
feed back device in coordination with reverse thrust to achieve the 
desired deceleration; 

3) To avoid using brake pressures that induce the skid; 

4) To decelerate mainly with reverse thrust using modulated braking 
for fine control and using nominal braking only when the maximum 
specified reverse thrust cannot decelerate the aircraft to the 
specified exit speed. 

This study was carried out for the TCV B737 aircraft and may be 
readily adapted to most conventional commercial aircraft. 



INTRODUCTION 


In order to reduce runway occupation time in rollout and turnoff it is 
desired to generate an automated uniform braking and reverse thrust guidance 
policy that will accommodate a large range of landing speeds, landing weights 
and runway conditions. To take full advantage of the aircraft capability in 
reverse thrust and braking it is necessary to base such a guidance policy on a 
fairly accurate model of the dynamic action of these control forces. 

No serious difficulty is encountered in generating a deceleration policy 
for dry, or even damp, runways in rollout and turnoff given the availability 
of brakes and reverse thrust (Ref. 1). However, for flooded runway conditions, 
hydroplaning in the ground speed range above 100 knots will negate the use of 
brakes. Moreover, the effective side force tire coefficient is rendered ineffec- 
tive on flooded runways above 80 knots. 

While it is possible to generate different guidance laws for diy and wet 
runways, it is preferable to devise a uniform policy, if at all possible. Such 
a policy is derived herein. 

Due to the stringent demands on time and space it is necessary to employ 
navigation aids such as MLS, buried magnetic leader cables, aircraft sensors 
and modem digital computer techniques to ensure safe and stable execution of 
the rollout, turnoff, and taxi commands. 

The detail description of the aircraft landing djmamics and navigation 
system (INS, MLS and the buried magnetic leader cable and aircraft sensors) 
is contained in Refs. 2 and 3, and is not repeated in this report. 

The report contains the derivation of the guidance laws and the results of 
the simulation for a variety of aircraft and runway landing conditions. 
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I. AIRCRAFT DECELERATION MODEL FOR BRAKING AND REVERSE THRUST 

For the purposes of modeling the deceleration of an aircraft in rollout, 
prior to turnoff, it is sufficiently accurate to treat the aircraft as a point mass 
in a single degree of freedom acting under the effects of drag, reverse thrust 
and braking. 

Let X be the position of the aircraft along the single degree-of -freedom 
coordinate, x its velocity and x ' , the acceleration. Then, we have for the 
deceleration of the point mass 


. . 1 .2 . . k 

X = - -T — p sC_ X - g - LL, H (x_) g + — 
2m ^ D ^R ^ m 


( 1 ) 


where 


R 


B 


H(iB> 


B 

k 

m 

g 

s 

P 


is the rolling friction coefficient 
is the tire braking friction coefficient 

is the heaviside operator initiated at the instant of braking 
command 

is the velocity at which braking is first permitted 
is the thrust 

is the aircraft point mass 
is the gravitational constant 
is the effective aircraft wing area 
is the mean air mass density 

is the drag coefficient in the landing configuration 


To integrate Eq. (1) we require the time history of the thrust, k (t), and 
the velocity x (t). 
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In the case of the thrust, we model the thrust system as a first order 
servo with a response time, r • Given a desired reverse thrust command. 


m 



-V 


-t/rx 

^ ) 


(2) 


where 


k is the thrust at the time of the reverse thrust command, 
o 

To integrate the drag term, we make use of the concept of averaging. 
The average squared velocity may be estimated from the assumption of 
a constant deceleration 



px (t) . 2 

\ X 

'Jx(o) 


d X 



(3) 


Thus we have 




^ t = ^ (t) + x(t) x(o) + x2(o^ (4) 

3 ^(t) - x(oJ) 


Employing Eq’s. (2) and (4) we have for the integration of Eq.(l) 


^ .2 . . .2 

' — T (x +X X +x )t-u„ gt 

6m ' o o ' ^ 


{ -Pb ® <*-V ) 

^ h C^^o" V T(l-e‘^/^)+kj^t) 


( 5 ) 






where t is the time of brake initiation. 
B 
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n, 


REVERSE THRUST AND BRAKING GUIDANCE FOR ROLLOUT 


AND TURNOFF 

Since the reverse thrust has a relatively long response time, it can 
not be used as a real time closed loop control variable. The brakes, how- 
ever, do possess good response characteristics. The philosophy used in 
generating a closed loop deceleration guidance law is to set the reverse 
thrust command at an average open loop deceleration level, well within its 
deceleration capability, and to utilize the brakes in a closed loop mode to 
realize a combined average constant desired deceleration. 

For this purpose, Eq. (5) may be used to solve for the desired reverse 
thrust, , that must be called for in order to reach the desired turn speed, 

, at the same time the aircraft reaches a desired position, . For this 
purpose, we assume that only the rolling friction is acting and that u = 0. 

To estimate the deceleration time, t , we use the assumption of constant 
deceleration, for which t is given by 


t = 


2 fx_ - X - - (t V /I.25TI 

L D o B ' s G D 


(6) 


+ X 

T o 


where D is a specified distance prior to turnoff when braking and reverse 
B 

thrust are not desired (in effect a safety margin) and t is a time over when 

s 

the brake pressure is increased to its nominal value. 


Finally, since t > > r » we have 


1 - 



1 


(7) 
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The desired reverse thrust, , is given by 

r 




‘'d == 


1-V. 


m (X “X ) - T k 
'To o 


psC 


D ,.2 -2 ^ ^ 

- (Xrr + X + X X ) + w 






( 8 ) 


where x^ is the desired turn speed 

X - X is the distance to go to the turnoff point 
D o 

W is the aircraft weight 

In the event that the called for reverse thrust exceeds the maximum reverse 

thrust, an estimate is made to see if the excess desired acceleration can be 

approximately achieved with additional braking. If it can, we set the command 

reverse thrust to and proceed. 

MAX 


The computation for the possible brake deceleration addition is given 
below: 




( 9 ) 


The possible brake friction coefficient is acceptable if; 
For dry and damp runways 


K, 


+ .02 s .4 


'^p W 

where . 02 is an arbitrary brake coefficient safety factor. 


(10a) 
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For wet and icy runways 


% w " :: — r 


(lOb) 


where we have that the hydroplane speed is given by 


X (KNOTS) = 9.V tire pressure 

H 


(Ref. 1) (11) 


In the event the test fails, we call for the next turnoff runway down the 
runway, and recompute the reverse thrust command utilizing the extended 
distance to go, x _ - x . 

O 


When the turning speed has been achieved, the thrust command is set 
to ensure a set taxi speed, x , opposing the deceleration forces of drag 

J. 

and rollii^ friction. Thus, we have for the taxi phase, the desired thrust 
for taxi, , is given by 


, 1 ^ • 2 
T 2 D TX ^ 


R 


W 


( 12 ) 


where the nominal value of 


^R 


.015. 


The braking law is designed to achieve and maintain a uniform average 
acceleration to reach the desired x^ , at the desired position, x^ , given 
the present acceleration velocity and position. The block diagram of the 
braking law is shown in Figure 1. 


We have for the average desired acceleration 



2 (DIST) 


(13) 
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where DIST is defined as 


DIST = 


^Distance to turnoff point - D 

a 

^7 meter (50) if DIST s 7 meters 


(13a) 


The braking feed back law is designed to maintain the constant desired 
deceleration. Let the error in the braking coefficient , , be given by 

X “ X 

0 < a = (13b) 

g 


The commanded braking coefficient , , is obtained from the solution of 

the differential equation 


■4- U = K^ ,K = 2 (14) 

dt '^e 

The commanded brake coefficient is limited to prevent excessive braking 
under dry conditions and to prevent skidding under wet conditions as follows; 

For dry runways 

Oiu (14a) 

L 
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For wet runways 

r 




0 for X > X 


H 


.014 X + 1, 
. 14 X + 2. 


^ .4 for X s X 


H 






where x is the hydroplaning speed. 
H 


(14b) 


The brake command is ramped to its nominal value during the first t^ 

seconds of brake initiation where t^ is a specified parameter ( t^ = 1 second). 

At t seconds after brake initiation, the control law is switched into the control 
s 

loop and the rate command integrator is initialized to the nominal brake coef- 
ficient, computed either in equation (lOa) or (10b). 
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'^nom 



^nom 

t 

s 

A 

B 


computed nominal braking coefficient 

time to apply nominal brake pressure 

limited brake pressure; for diy runway set equal to . 4 

for wet runway set equal to . 0 14 x + i 

. 14 X + 2 

if wet runway and ground speed greater than hydroplane speed 
initialize integrator at t^ to 


Figure 1. Closed Loop Braking Law 










m. STEE RING GTHDANCE FOR ROLLOUT TURNOFF AND TAXI 


The object of the steering guidance logic is to direct the aircraft to 
follow a given sequence of straight lines and curves from touchdown to turn- 
off to taxi in a short time span so as to safely reduce runway occupancy. In 
a manner completely equivalent to the in flight way point RNAV equations , 
we fix a set of runway coordinate points consisting of the beginning and end of 
each segment. See Figures 2a and 2b. 

The guidance parameters to be determined on each segment consist of 


the following: 


1) 

Cross track error 

(CTE) 

2) 

Cross track rate error 

(CTRE) 

3) 

Track angle error 

(TAE) 

4) 

Yaw rate error 

(YRE) 

5) 

Distance to go to the next 
way point 

(DTGO) 

6) 

Time to go to the next 
segment 

(TTGO) 
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Figure 2 (b). Runway (Inertial) Coordinates and Aircraft Position, Velocity Vectors (R,R) 
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From the navigation system we may obtain the best estimate of aircraft 

A A A A A 

position,velocity and acceleration in runway coordinates (x , y , x , y , x , 

A A A 

y ) as well as the aircraft yaw and yaw rate ( j/) , ) • The on board com- 

puter contains the coordinates of the way points ^iVX(I), WY(I)^ , the yaw 
heading angle ( W0(I) ) at the ith segment, and the length of each segment, 
LSG(I). 

The guidance parameters are detailed below, (1) for straight line seg- 
ments, (2) for circular segments and finally (3) for centered parabolic seg- 
ments. (See Appendix A). 

The aircraft deviations from the desired location on runway track (TX , 
TY) is not necessarily computed with respect to the eg of the aircraft, but 
rather with respect to a sensor on the aircraft center line, (AXP , AZP), 
whose location is chosen to be in a logical belly location of the aircraft. 

(See Figure 3). 

The coordinates of the aircraft sensor in runway coordinates is given by 


XRP = X + cos ill)) AXP 

A A 

YRP = y + sin ( j/) ) AXP 


(15) 


The coordinates of sensor velocity in runway coordinates is given by 


XRP 

YRP 


X - sin ( 0 ) AXP 

A A 

y + cos (ij) ) AXP 


r 

A 

r 


(16) 
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Figure 3. Landing Gear Geometry 
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(1) straight Line Segments 


(17a) 
(17b) 
(17c) 
(17d) 

(17e) 
(17f) 

where 

LSG(I) = lei^h of the Ith straight segment 
RAD = 180 . /tt 

WX(I) , WY(I) , Wi/)(I) are defined at start of each segment. 

The time to go equation is common to all segment types. 

A test is made to switch to the guidance parameter logic for the next 
segment whenever 


GTE = (^YRP - WYa)^cos (^W0(I)^ - (xRP - WX(I)^ sin (w 0(1)^ 
CTRE= YRP cos (w 0(1)^ -XRP sin (w w) 

TAE = 


YRE = r 


180 

7T 


DTGO = LSG(I) - (yrv - WY(I)^ sin (w 0(1)) - (xRP - WXa)) cos (w 0(1)) 


TGO = DTGO /( ) 


1/2 


TGO ^ .05 second ^18) 

(2) Circular Segments 

The desired runway track coordinates for the circular turn are computed 
from the intersection of the circular track and the radial vector from the center 
of the turn circle through the best estimate of the aircraft position. 
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The coordinates of turn circle center are given by 

CTX(I) = WX(I) - SGN(I) RT(I) sin (w ^(I) ) 
CTY(I) = WY(I) + SGN(I) RT(I) cos (w j/)(I) ^ 

when SGN(I) = +1. for a right hand turn 

= -1. for a left hand turn 


and RT(I) = radius of the turn 

The coordinates of the radius vector from the turn center to the aircraft 
sensor are given by 


XS = XRP - CTX(I) 

YS = YRP - CTY(I) 

GTE = SGN(I) [ RT(I) - + YS^) ] 

• • I ^ 2 ^ 

CTRE= -SGN(I) (XS XRP+ YS YRP) / (XS + YS ) 



TAE 


YRE 


TAN 


-1 


XS cos (0) + YS sin (0) I 180 


YS cos (j/j) - XS sin (ij)) 

[ r - (x^ + y^)^ SGN(I)/RT(I) ] ^ 


( 20 ) 

(21a) 

(21b) 

(21c) 

(21d) 


DTGO= RT(I) ij)^ - W ?/j(I+l) 


(21e) 


where 


- W ib(l+l) = tan 
Lto 


-1 ) - sin (W 0(1+1) YS - cos (W ^(I+l) XS 
-cos (Wij){I+l) YS + sin (W j/j(I+l) XS 


(21f) 
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The concept of a parabolic turn segment is introduced in order to 
attempt to reduce the severity of the initiation of the yaw rate onset during 
a fast turnoff. It represents an "easy on" and "easy off" for a turn and 
reserves the maximum turn rate for the center of the turn. The theory and 
derivation is contained in Appendix A. The guidance parameter equations are 
detailed below. 

(3) Centered Parabolic Segments 

The desired position on the parabolic track is most conveniently de- 
scribed in a cartesian coordinate system whose origin is at the middle of the 
parabolic turn as described and illustrated in Appendix A. The desired position 
on the parabola requires the solution of a cubic polynomial. The equations must 
be solved in real time and are outlined below: 

The vector from the start of the parabola to the aircraft sensor in runway 
coordinates is given by 


XS = XRP - WX(I)^ 

YS = YRP - WY(I) ) 


( 22 ) 


Let 


a 


w ^(i) + w 4){i+i) 

2 




W - W 

2 


RP(I) = 


cos 0 ! ^WX(I+1) - WX(D) + sin (WY(I+1) - WY(p') 
2 tan 


(23) 


XPAI 


-SGN(I) RP(I) tan ^ 
1/2 XPAI tan $ 
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YPAI 



The coordinates of the aircraft sensor in the centered parabolic system are 

XPA = XPAI + XS cos Of + YS sin o 
YPA = YPAI - XS sin a + YS cos a 

The desired position on the parabolic track in centered coordinates is given by 


where: 


XD 

AQ 

BQ 

YD = 
e: 

1/2 

RP(I) 

SGN(I) 

AQ = 

-HQ + SQ 


BQ = 

-HQ - SQ 


SQ = 





HQ = 

1/2 XAQ 


DISC = 

4 XAP + 27 XAQ 

XAQ = 

-2 RP^(I) XPA 

XAP = 

2 RP(I) rRP(I) - YPA 


The desired yaw heading at the desired track position in the parabolic system 
is given by 


sin i}) 


DP 


xr> 

(XD + RP ) 


, _ RP 

^DP - — ^ 

(KB + RP ) 
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The desired heading in the runway system is given by 


D 


+ OL 

^DP 


The desired yaw rate on the track is given by 


D 


/'?2 

ix +y j SGN(I) 

T“ 

2 2 ^ 
(XD + RP ) 


In terms of the above we may now compute the required guidance 
parameters 


GTE = (YPA - YD) COS ( - (XPA - XD sin 


CTRE = YRP cos 0^ - XRP sin 0^ 


-rD(YPA sin + XPA cos ) 


-r (- 1/2 SGN(I) RP(I) sin (1 + ^ ) ) 

COS 0 


DP 


TAE = tan 


-1 


(sin 0 cos 0p - cos 0 sin 0^) 


(cos 0 cos 0j^ + sin 0 sin 0^^) 


180 

TT 


YRE = r - r 


D 


sin 0 


DTGO = 1/2 LGS(I) - SGN(I) RP(I) ( + In 


1+sin 0 


2 , 

cos 0 


DP 


cos 0 


S ,) 


DP 
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(27) 

(28) 

(29a) 

(29b) 

(29c) 

(29d) 

(29e) 
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with the buried magnetic leader cable in operation the cross track 
error and the track angle error may be taken directly from the voltage out- 
put of the three magnetic coils situated at sensor locations. The details of 
theoiy of the magnetic cable and sensor are described in Refs. (2) and (3). 

Let the voltages from the output of the three coils be v , v , and v , 

X ’ y" z 

then we may obtain directly 

V 

CTE = (ZP-Z^j^) ^ (30) 

y 

V 

TAE = “ (31) 

V 

y 

when ZP - is a known height of the pickup above the buried cable. 

WIRE 

The steering law is a slight modification of the Boeing B737 rollout 
equation which uses a nose wheel command geared to the rudder signal by 
the ratio of 


NW 


= - 7/26 6 


COM 


RUD 


COM 


The block diagram for the steering law is given in Fig. 4. 


(32) 
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Figure 4. Rollout and 
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B7 = 2. 


Steering Command 





IV. DESCRIPTION OF THE NAVIGATION SYSTEM 


The navigation system of the B737 consists of a stable platform which 
provides the rigid body inertial altitude and body angular rates. These are 
noisy and are subject to drift. Body mounted accelerometers provide a con- 
tinuous measurement of the accelerations due to the external forces. These 
measurements are noisy and contain errors due to misalignment and scale 
factor errors. 

Ref. (2) describes both a Kalman type filter and a constant gain 
complementary filter, either of which may be used to filter out the noisy 
measurements and provide for corrections so that the body mounted accel- 
erometers may be integrated in an inertial reference frame to provide position, 
velocity acceleration, aircraft altitude and angular rates to the RNAV and 
displays. Fig. V is a block diagram of the Rollout and Turnoff Navigation and 
Guidance System. 

Upon touchdown the major error source in the aircraft position is the 
down range position error due to the MLS range bias which cannot be re- 
moved by filtering. The range bias produces errors in the guidance signals 
which results in large tracking errors during the turn. This is illustrated 
in the first simulation run in Fig. VI. To eliminate this error source we 
recommend that the aircraft should receive a short radio blip signal as it 
passes a known position along the runway during rollout. By correcting the 
X coordinate to the known runway position, the offset problem disappears and 
good turnoffs may be realized. 

The update equations are as follows : 

Let the distance from the second way point to the blip location be 
BLIPX. Then 

XDIST = LSG(l) - BLIPX 
X (t ) = XDIST 

^ blip' 


(33a) 

(33b) 
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The use of a buried magnetic leader cable as described in Refs. (2) and 
(3) is imder test at Langley Field and Wallops Island. Some results from these 
tests are given by Ref. 4. These tests indicate that while the use of sensors 
may provide a usable signal for the cross track error, the possibility of ob- 
taining the track angle error (TAE) directly is in doubt. In any case, with the 
major source of error in the MLS measurement removed by the preceding 
BLIP technique, it should be possible to utilize MLS data to estimate TAE. 
Additional tests and simulation are recommended. 
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Figure 5. Rollout and Turnoff Navigation and Guidance System 






V. SIMULATION RESULTS 

This section contains the computer results of the simulation sensitivity 
study of the Langley TCV B737 aircraft in rollout, turnoff, and taxi. The 
runway geometry simulated was the Wallops Flight Center high speed exit 
on runway 22. 

The factors affecting the study that were systematically varied are as 
follows : 

1) The effect of MLS range bias error on lateral ’'stand off" in turns 

2) The effect of aircraft speed in turns 

3) The effect of winds, varying both the direction and magnitude 

4) The effect of diy vs wet runways 

5) The effect of magnetic sensor noise and bias 

6) The effect of filtering methods on performance 

7) The effect of turn radius 

8) The effect of first exit distance on runway occupancy time 

9) The effect of landing weights 

10) The effect of parabolic vs circular turns on lateral acceleration 
in turns. 

In order to illustrate the results of the study CALCOMP plots of the 
computer runs are shown in Figures 6 through 23. Table I contains a list 
of the parameters used in the simulation study for each plot. 

Each simulation run consists of 12 plots of significant variables vs time. 
The items plotted are: 

1) Thrust magnitude (Newtons) 

2) Percent brake pressure command 

3) Along track deceleration (g's) 

4) Ground speed (m/sec) 
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TABLE I 

ROLLOUT SENSITIVITY RUNS 


RUNWAY TURN TAXI 

COND. SPEED SPEED 

KN KN 


DIST. TO I WINDS 
FIRST I speed/ 
/der 



KALMAN 

COMP 

NOISE 

AND 

BIAS 


-j 




DRY 

65 

60 

11 

65 

60 

M 

65 

60 

M 

45 

40 

11 

65 

60 

11 

75 

70 

M 

45 

40 

11 

75 

70 

11 

75 

70 

tf 

65 

60 

WET 

65 

60 

DRY 

65 

60 

11 

65 

60 

M 

65 

60 

H 

65 

60 

11 

65 

60 

DRY 

65 

60 


10/45 


10/45 


548.6 

548.6 


548.6 

914.4 

548.6 


548.6 

548.6 

548.6 


40823. 3 
40823.3 

ii 


40823.3 

31751.5 

40823.3 


Excessive 

Noise 

YES 

















5) 


Yaw (Degrees) 

6) Cross track error (meters) 

7) Nose wheel steering command (degrees) 

8) Error in estimates of runway position coordinates for the 
navigation filter used in simulation run 

9) Error in estimates of velocity components in the runway 
coordinate frame for the navigation filter used in the 
simulation run 

10) Track angle error (degrees) 

11) Cross track velocity error (m/sec) 

12) Cross track acceleration (g’s) 

A brief discussion of the significant features of Figures 6 through 
21 follow. 

Figure 6 contains a t5^ical automated rollout and turnoff when no 
information is available to the navigation system to eliminate the MLS range 
bias in the DME signal. The rollout is completed in an acceptable manner, 
however, the cross track error (Plot #6, Figure 6) shows a fixed "stand off" 
resultir^ from the down track error in position (Plot #8, Figure 6) settling 
into a cross track error in the turn. The fix for this effect is to give the 
navigation system a more accurate along track position measurement early 
in the rollout so that the system knows where it is relative to the beginning 
of the turn. Figure 7 and all subsequent figures contain this blip information 
and a readjustment can be noted in x-position error (Plot #8, Figs. 6 - 21) 
shortly after touchdown. It may be noted that the large "stand off' error in 
cross track is no longer present in these figures. 

Figures 7 and 8 illustrate the performance of both the Kalman and 
complementary filters as navigation aids when no noise or bias is present in 
the magnetic leader cable. Response to cable noise and bias due to conducting 
reenforcing rods in the runway pavement construction is mainly present in the 
nose wheel steering command channel (Plot #7, Figs. 9 - 19). In contrast to 
Figs. 7 and 8 one may examine Fig. 20 in which excessive noise drives a 



nose wheel command chatter. The aircraft is moving too rapidly to follow 
the nose wheel motion, however the effect is to produce an unacceptable 
nose wheel response. The results of the tests in Ref. 4 would indicate that 
the yaw deviation signal coming from the magnetic cable sensor may have to 
be eliminated and possibly replaced by information from the INS and pre- 
stored information regarding the runway center heading. Additional study 
will be required to solve this problem. 

Figures 9, 10, 11, 12, and 13 illustrate the performance of the auto- 
mated rollout and turnoff as a function of turn speed using either Kalman or 
the complementary filters. The main features to be noted are that the cross 
track acceleration (Plot 12, Figs. 11 and 13) for a 75 knot turn produces a 
maximum lateral acceleration of 0.3 g's which borders on the uncomfortable. 
This is due to an S shape reversal in the turnoff pattern (Figure 2). An 
acceptable sustained maximum lateral acceleration is on the order of 0.2 g’s 
during a steady turn. Figure 14 illustrates the effect of attempting to alleviate 
the severity of the turn rate onset using a "centered” parabolic turn. It is 
evident that no significant improvement is forthcoming from this technique. 

The improvement was small because the "centered" parabola path did not 
differ significantly from the circular path for this 30^ turn. It is possible 
that by "easing in" the turn signal in the guidance law loop may help to 
eliminate this effect. Additional work remains to be done in this area if 
75 knot turns are to be mandated for runway occupancy reduction. 

It is significant to note that turnoffs at speeds of 45 and 65 knots may 
be executed in this automated mode. It is questionable whether a manual 
operated system could perform as well at these speeds. 

Figure 22 illustrates the effect of turnoff distance on runway occupancy. 
Contrasting Fig. 10 and 22 both of which call for the same turnoff speed at 
different turnoff distances (Figure 22 turnoff distance is 1300 feet longer) show 
almost identical behavior except that there is practically no braking for the 
longer turnoff distance since the reverse thrust is able to achieve the turn speed 
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at the proper turn distance without additional braking. The runway occupancy 
increases from 15 seconds for the shorter turnoff runway as compared with 
20 seconds for the longer turnoff runway in Fig. 22. 

The effect of wet runways is illustrated in Fig. 16. This case is to 
be compared with Fig. 9 both of which call for a turn speed of 45 knots. The 
wet runway case cannot achieve the deceleration in time to make the turn at 
first turn opportunity since the ground speed exceeds the hydroplane speed 
and no braking is permitted. By calling for the next turn opportunity, the 
wet runway case is able to achieve the turn speed at just under 25 seconds and 
is able to execute a safe exit at 45 knots. 

Figure 18 illustrates the effect of turn radius on the cross track accel- 
eration. By contrasting Fig. 17 with Figure 10 we obtain the expected results. 

The sustained lateral acceleration in the steady 65 knot turn at a radius of 548.6m 
is reduced from . 2 g’s to . 11 g’s at 65 knots and 914. 4m radius. 

The effect of winds in varying directions is illustrated in Figures 18 and 
19. These runs are to be contrasted with Fig. 10. All three call for a turn 
speed of 65 knots. The effect of winds is to increase the maximum cross track 
error from 4 meters for the case of no winds to 6 meters for 10. knot winds. 

The cause of the overshoot is that winds provide for additional yaw rate as the 
aircraft heading changes making it more difficult for the rollout control law 
to maintain the desired headir^ in a turn. It is possible that feeding in infor- 
mation on wind direction could provide for a better control in cross track error 
during the turn in a wind condition. Additional study is required in this area. 

Figure 21 shows the effect of aircraft weight. The major effect is to 
reduce the braking required for deceleration, since the reduced weight (40823. 3 kg 
in Fig. 10 is reduced to 3175L5kg in Fig. 21) permits the reverse thrust to 
achieve the turn speed of 65 knots without much brake assist. 
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VI. CONCLUSIONS AND RECOMMEISTDATIONS 


In general, it is concluded that runway occupation time can be 
significantly reduced by use of automated rollout and turnoff capability. 

The aircraft can be clear of the runway within 20 seconds for diy runways 
and 30 seconds for wet runways. It is further concluded that turns of 
33. 4 m/s (54 knots) at 548. 8 in. (1800 ft) radii or laiger can be automatically 
and safely achieved. 

It is recommended that rollout and turnoff guidance laws to reduce 
cross track errors due to winds be investigated. 

It is also recommended that "slow in" logic be investigated to ease 
the onset of turn acceleration. 
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Turn Speed ■ 45 knots 
Taxi Speed ® 40 knots 
No Winds 

Turn Radius » 548 . 6 m 
Landing Wt. ■ 40823. 3 kg 
Kalman Filter Used 
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Figure 13. Rollout and Turnoff on Dry Runway - Noise and Bias Acting 
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Figure 14. Rollout and Turnoff on Dry Runway - Noise and Bias Acting 
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Figure 16. Rollout and Turnoff on Wet Runway - Noise and Bias Acting 
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Figure 18. Rollout and Turnoff on Dry Runway - Noise and Bias Actir^ 

Turn Speed * 65 knots 

Taxi Speed 60 knots 

Winds » 10 knots at 45° heading 
Turn Radius » 548.6 m 

Landing Wt. » 40823. 3 kg 

Kalman Filter Used 












58 




Turn Speed * 65 knots 

Taxi Speed » 60 knots 

Winds « 15 knots at 0 heading 
Turn Radius * 548.6 m 

Landing Wt. » 40823. 3 kg 

Kalman Filter Used 
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Figure 22. Rollout and Turnoff on Dry Runway - Noise and Bias Acting 
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APPENDIX A 


RNAV EQUATIONS FOR PARABOLIC TURNS 
IN ROLLOUT AND TURNOFF 



RNAV EQUATIONS FOR THE 
GENERAL RUNWAY CURVE 


Let the runway coordinate system have its origin at the apron of the 
landing strip, with the x axis pointing north, the y axis due west, and z axis 
down. Yaw is measured clockwise from the north pointing x axis. 

Let the prescribed curve, describing the required turnoff path, be given 
in the most general form 

N 

y = f (X ) = y a. x^ 

■^D ' d' a i D 

i=0 

The required yaw for any point on the curve is given by 


dy 


N 


tan0^ = 


D 


dx 


D 


= y i a. x^ 
L i D 

i=0 


i-1 


(A2) 


To find the lateral deviation from the prescribed curve at some x, y 
point, not on the curve, we have 


6 = (y - Yq) cos - (X - Xj^) sin 0^ 


(A3) 


It is necessary to determine x^ , y^^ , and , given x, y and the 
coefficients, a. , defining the curve. 

We have 


X - X 

y - y 


D 

D 


1 

'^D/iyD 


N 


-I 


i i-1 
1 a Xj^ 


i=0 


(A4) 
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D 


Substituting Eq^Al) for 
that must be solved, 


D 


into Eq.(A4) we have for the polynomial in 


N 




- (y 


-I 




N 
0 


i-l 


(AS) 


It is plain that for a given curve of order N , the polynomial to be solved 

for in X is of order 2N - 1. For N > 2 this will require some Newton 
D 

iteration for real time processing of the RNAV equations. Such a computing 
burden is to be avoided, if possible. 

For N = 2, the most general conic section requires a cubic. For the 
special case of the circle, only a quadratic equation need be solved; which 
accounts for its use in inflight RNAV. However, if we desire to avoid exces- 
sive abrupt lateral accelerations we must ease into the turn and other conic 
sections should be examined. 

For the centered parabola we obtain a cubic which we may guarantee has 
only one real root for reasonable deviations from the prescribed curve. 

We define a centered parabola as 

"d ' il" (A6) 

o 

where R is the minimum radius of curvature at the point ( o , o ) , and 
o 

sgn (0^-0 ) is the sign of the difference between the required yaw at the end 
of the turn and required yaw and the beginning of the turns. 

The angle between the tangent to the centered parabola and the v axis is 
given by 


du 

Tv 


D 

D 


= tan /3 


D 


sgn ( 0j - 0^) Vj^ 


R 


(A7) 
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To obtain the lateral deviation from the prescribed curve at a 

point (u, V), not on the curve, we have 


+ 2R^ (r^ - sgn ( - 0^) u ) Vj^ - 2R^^ v = 0 (A8) 

The solution of Eq. (iA6) will have one real root provided 

2 

sgn ( ) u < “ ( 27/4 R^ v ) (A9) 

For reasonable distances from the prescribed curve, only one real solution 
exists and it is given by 

= a/ I A i + y 1 B (AlO) 

where 

3 2 

3 2 

and 


2 

Q = - 2 R V 
^ o 

P = 2 R^ - sgn ( - l/j^) u) 


(AlOb) 


The corresponding value of u^^ is given by Eq. (A6). 
Uj, = sgn (((If- 1(1^) vP/ 2 H^ 


(A6) 
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To ejqpress u , v in terms of the runway coordinates, we have 



r '‘’o'"'’! 'i 

1 I 

2/0 A 

itan ( — 5 — ) 


where 


X , y are the initial runway coordinates at the start of the turn, and 
o o 

QL is the average yaw angle. 


To obtain the cross track error we have from Eq. (3) 
CTRE =6 = (11 - Uq) cos 0^- {V - sin 


where 


cos = 


R 

o 

2 TT 


singj, = 


sgn 

2 2 i 

<Ro 


u, V are obtained from the runway coordinates ( x, y ) by Eq. (11) and u^^ , 

V . cos it t and sin ib are the solutions of the cubic given by Eqs. (AlO), 
(A 6 ), and (A;l2a). 

The track angle error is obtained fix>m 
TAGE = ^ - 


(All) 


(Alla) 


(A12) 


(A12a) 


(A13) 


70 



where 0 is the true yaw 


and 




D 


+ a 


A useful formula for the track angle error is given by 


sin 0 cos 0 - cos 0 sin 0 

TAGE = tan"^ ( ; : ^ ; : ) 

\ cos 0 cos 0j^ + sin 0 sin 0^ / 


The required yaw rate is given by 

//' 2 2 2 N 

*'’d ° ''g < '<’f ■ ‘(’o >/0\ ^ ''d* ''%) 

2 2 2 

The expression in the denominator, ( + v^ ) ), is the radius of cur- 

vature of a point on the parabola ( u^^ , v^^). The sharpest turn is seen to lie 
in the center of the turn as desired. 

The required cross track error rate is given by 


6 = cos y - V cos 8^^ - sin 0^ x - u sin 3^^ 0^ 

• ^o 2 

+ 0D — sgn ( 0^ - 0^) sin 0^ (1 + sec 3^^) 


To obtain the distance to go to the end of the turn we have 




tan 3j sec 0^ + In (tan 0^ + sec 0^} 


DSTGO = i sgn ( 0^ - 0^) R^ / 


-tan 3j^ sec 0^^- In (tan 3^+ sec 3^) 

V. 


(A13a) 


(A14) 


(A15) 


(A16) 


) (A17) 
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where 


It 


2 

r 

sgn (^j - il)^) 

tan 0^ 

R 

o 

1 

R 

[sec 0^ = 

o 

2 2 4 


(^0 > 


(A 17a) 


It is convenient to relate the maximum radius of curvature, R^, to the 
equivalent circular radius of curvature between the same two end points. 


We have for the parabola 

cos a (Xj " x^) + sin a (y^ " y^) 


R 


2tan( —2 ) 


(A18) 


For the circle through the same end points 

(cos - cos j/)^) (Xj - x^) - (sin - sin (y^ - y^) 


R 


2( cos - cos 0^) (sin j/)^ - sin tlj^) 


(A19) 


Figure 23 illustrates the two curves through the same end points. 


72 




REFERENCES 


1. Horne, Walter B. ; McCarty, John L. ; and Tanner, John A. ; "Some 
Effects of Adverse Weather Conditions on Performance of Airplane 
Antiskid Braking Systems". NASA TN D-8202, 1976. 

2. Pines, S. ; Schmidt, S. F. ; and Mann, F. ; "Automated Landing, Rollout 
and Turnoff Using MLS and Magnetic, Cable Sensors". NASA CR-2907, 
1977. 

3. Olsen, K. ; Ventola, R. ; and Pines, S. ; "Tests of Magnetic Leader Cable 
and Sensors as a Rollout Navigational Aid". Analytical Mechanics 
Associates, Inc., Report No. 77-22, Nov. 1977. 

4. Hueschen, R. M. ; Creedon, J. F. ; Bundick, W. T. ; and Yovmg, J. C.; 
"Guidance and Control System Research for Improved Terminal Area 
Operations, " Presented at the 1980 Aircraft Safety and Operating 
Problems Conference at Hampton, Virginia, November 5-7, 1980. 



1. Report No. 

NASA CR-3451 


2. Government Accession No. 


3. Recipient's Catalog No. 


4 Title and Subtitle 

TERMINAL AREA AUTOMATIC NAVIGATION, GUIDANCE, AND CONTROL 
RESEARCH USING THE MICROWAVE LANDING SYSTEM (MLS) 

PART 1 - AUTOMATIC ROLLOUT, TURNOFF, AND TAXI 


5. Repoa Date 
August 1981 


6. Performing Organization Code 


7. Author{s) 

Samuel Pines 


8. Performing Organization Report No. 

AMA 80-1 


to. Work Unit No. 


9. Performing Organization Name and Address 

Analytical Mechanics Associates, Inc. 
17 Research Road 
Hampton, Virginia 23666 


11. Contract or Grant No, 

NASl-15116 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


13. Type of Report and Period Covered 
Contractor Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Langley Technical Monitor: Richard M. Hueschen 
Final Report 


16. Abstract 

This report contains the results of a study developed for the Langley TCV B-737 
program to apply existing navigation aids plus magnetic leader cable signals and develop 
braking and reverse tiarust guidance laws to provide for rapid automated rollout, turnoff, 
and taxi to reduce runway occupation time for a wide variety of landing conditions for 
conventional commercial-type aircraft. 

Closed loop guidance laws for braking and reverse thrust are derived for rollout, 
turnoff, and taxi, as functions of the landing speed, the desired taxi speed and the distance 
to go. Brake limitations for wet runway conditions and reverse thrust limitations are 
taken into account to provide decision rules to avoid tire skid and to choose an alternate 
turnoff point, farther down the runway, to accommodate extreme landing conditions. 


17. Key Words (Suggested by Author(s)) 
Rollout 
Turnoff 
Braking 
Reverse thrust 
Steering guidance 


18. Distribution Statement 


Unclassified - Unlimited 


Subject Category 08 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price 

Unclassified 

Unclassified 

79 

A05 


For sale by the National Technical Information Service, Springfield. Virginia 22161 


NASA-Langley, 1981 




